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The low-temperature properties of CeOBiS 2 single crystals are studied by electrical resis¬ 
tivity, magnetization, and specific heat measurements. Ce 4/-electrons are found to be in a 
well-localized state split by crystalline-electric-field (CEF) effects. The CEF ground state is 
a pure J^ = +1/2 doublet, and excited doublets are located far above. At low temperatures 
in zero field, we observe pronounced - log T divergence in the specific heat, revealing the 
presence of quantum critical fluctuations of 4/ magnetic moments near a quantum critical 
point (QCP). Considering that CeOBiS 2 is a nonmetal, this phenomenon cannot be attributed 
to the competition between Kondo and the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter¬ 
actions as in numerous /-electron-based strongly correlated metals, indicating an unconven¬ 
tional mechanism. We suggest that CeOBiS 2 is the first material found to be located at a QCP 
among geometrically frustrated nonmetallic magnets. 


The - log T divergence of specific heat has been observed in numerous /-electron-based 
strongly correlated electron systems.The - log T behavior is considered to be a hallmark 
of non-Fermi-liquid (NFL) behavior realized around a quantum critical point (QCP) or a 
r = 0 phase transition, where quantum fluctuations of magnetic moments dominate. Since 
the realization of such a QCP requires the balance between Kondo and the Ruderman-Kittel- 
Kasuya-Yosida (RKKY) interactions, fine tuning of the relevant control parameter by chemi¬ 
cal doping is needed in many cases. Since chemical doping inevitably introduces randomness 
into a crystal, a QCP realized in a pure compound is beneficial for research to clarify the 
mechanism. However, the numbers of such examples (including CeNi 2 Ge 2 ^ and jS-YbAlB 4 '^) 
are limited. In this letter, we report the finding of - log T divergence of the specific heat in 
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nondoped CeOBiS 2 . Since this material is a nonmetal, the Kondo effect (a Fermi surfaee ef- 
feet) cannot be realized in principle. Therefore, an unconventional mechanism is necessary 
to aeeount for this behavior. 

CeOBiS 2 is a mother phase of a reeently found BiS 2 -based superconductor.^^^ As shown in 
Fig. 1, the crystal structure consists of an alternating staeking of BiS 2 and CeO layers along 
the c-direetion. In CeOi_j:F;tBiS 2 , the partial replacement of O by F provides eleetron doping 
into the BiS 2 layers and induces superconductivity with the transition temperature of 2.3 
K (under applied pressure, it increases to as high as 6.7 K^°). In view of the two-dimensional 
(2D) nature of the crystal structure, it is intriguing to compare the physieal properties of 
CeOi_;i.Fj.BiS 2 with those of euprates'^ and Fe-based superconductors.^^’Single-crystal 
studies have shown that CeOi_j;F;cBiS 2 has strongly anisotropie supereondueting properties. 

In addition, the possible eoexistenee of superconductivity and Ce 4/-electron ferromagnetism 
in this system has been pointed out reeently using polycrystalline samples.*^’ These findings 
have motivated us to investigate in detail the 4/-eleetron magnetism using single crystals. In 
this study, we investigate the physical properties of undoped CeOBiS 2 single erystals. This 
is the first report on the magnetism of 4/ eleetrons in the series of LnOi_j:F;tBiS 2 {Ln\ rare 
earth) supereonduetors using single crystals. 

Single crystals of CeOBiS 2 were prepared by the CsCl flux method using Ce chips (99.9%) 
and powders of Ce02 (99.99%), Bi 2 S 3 (99.999%), S (99.9%), and CsCl (99.9%). The excess 
CsCl flux was removed using H 2 O. The single crystals are plate-shaped parallel to the c-plane 
with typieal dimensions of 1 x 1 x 0.2 mm^ (see inset of Fig. 2). The powder X-ray diffraetion 
speetra shown in Fig. 2 confirm the LaOBiS 2 -type strueture with PAInmm symmetry. For the 
determination of detailed structural parameters, single-crystal X-ray analysis was performed 
using a Rigaku Mercury diffractometer with graphite monochromated Mo-Ka radiation. A 
seleeted small single erystal with dimensions of about 0.10 x 0.10 x 0.05 mm^ was mounted 
on a glass fiber with epoxy. Structural parameters refined (with the occupation ratio for each 
site fixed to 1) using the program SHELX-97’^ are summarized in Table I. The lattiee pa¬ 
rameters agree well with the reported values.*^ DC magnetization (M) measurements were 
earried out in a Magnetic Property Measurement System (MPMS; Quantum Design (QD)) 
down to 2 K and up to 7 T. Speeific heat (C) measurements were performed using a quasi- 
adiabatic heat-pulse method with a Physieal Property Measurement System (PPMS; QD) and 
a dilution refrigerator down to 0.2 K up to 8 T. Considering that the measured single crys¬ 
tal includes nonmagnetic Ce02 and CsCl phases (see Fig. 2), a sample mass correction is 
made by multiplying the measured values of M and C by a eorreetion faetor of p = 1.12. p 
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Fig. 1. (Color online) Tetragonal crystal structure {P4lnmm) of CeOBiSi, consisting of an alternating stack¬ 
ing of BiS 2 and CeO layers (depicted using VESTA®). A unit cell is indicated by solid lines. 


is determined by erystalline-electric-field (CEF) model fitting to the M(T) data (see below). 
Eleetrieal resistivity was measured by a standard AC four-probe teehnique. 

In Fig. 3, the temperature dependenees of the in-plane resistivity Pa are shown for sin¬ 
gle erystals of CeOBiS 2 , EaOBiS 2 , and EaOo. 5 Fo. 5 BiS 2 . The Pa(T) data of CeOBiS 2 and 
EaOBiS 2 show semiconducting behaviors and are three orders of magnitude larger than 
those of electron-doped EaOo. 5 Fo. 5 BiS 2 . Such extremely weak conduction of CeOBiS 2 and 
EaOBiS 2 attests to the low impurity (or defect) concentration and high quality of the sin¬ 
gle crystals. If the semiconducting behavior of CeOBiS 2 is analyzed tentatively using the 
activation-type relation pa = poexpiAIlk^T) [see Fig. 3 (b)], the energy gaps Aiow = 0.93 K 
below 5 K and Ahigh = 95 K above 100 K are obtained; see Ref. 18 for a similar analysis of 
the semiconducting behavior observed in FaOo. 5 Fo. 5 BiS 2 polycrystals. Since these values are 
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Fig. 2. (Color online) X-ray powder diffraction pattern for crushed single crystals of CeOBiS 2 using a Cu- 
Kal radiation. The Miller indices of the tetragonal symmetry P4lnmm for main peaks are shown. There are a 
few peaks attributed to the remaining Ce02. (inset) Picture of CeOBiS 2 single crystal. 


Table I. Atomic coordinates and thermal parameters of CeOBiS 2 at room temperature determined by single¬ 
crystal X-ray measurements. R and wR are reliability factors and Bgq is the equivalent isotropic atomic displace¬ 
ment parameter. Standard deviations in the positions of the least significant digits are given in parentheses. 


P4lnmm (11129) a = 3.984(6) A, c = 13.49(2) A, T = 214.2(6) A^ 
(origin choice 2) Position 


Atom 

site 


X 

y 

z 

Beq (A2) 

Ce 

2c 

{4mm) 

1/4 

1/4 

0.0919(2) 

0.35(7) 

Bi 

2c 

{4mm) 

1/4 

1/4 

0.6272(2) 

0.76(7) 

S(l) 

2c 

{4mm) 

1/4 

1/4 

0.3817(11) 

1.2(3) 

S(2) 

2c 

{4mm) 

1/4 

1/4 

0.8126(9) 

0.9(3) 

O 

2b 

(4m2) 

3/4 

1/4 

0 

0.4(6) 

R = 8.39 %, wR = 

13.75 % 






much smaller than the band gap energy of 0.4 eV obtained from a band strueture ealeulation 
on LaOBiS 2 ,'^ the observed Aiow and Ahigh may be attributed to minor impurity (or defeet) 
eonduetion. 

For the eleetron-doped LaOo.sFo.sBiSi, Pa(T) shows a supereondueting transition at Tc ^ 4 
K. When T > T^, Pa{T) shows weak but elear metallie behavior, in eontrast to the semieon- 
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Fig. 3. (Color online) (a) Temperature dependence of the in-plane resistivity pa for single crystals of L«OBiS 2 
{Ln - La, Ce) and LaOo. 5 Fo. 5 BiS 2 at 0 T. (b) 1/T dependence of logpa for CeOBiS 2 . The solid lines are simple 
activation-type fitting curves in the low- and high-temperature regions. 


ducting behavior observed in polycrystals. 

Figure 4(a) shows the temperature dependences of the DC susceptibility of CeOBiS 2 mea¬ 
sured in 0.01 and 7 T for // // a. No anomaly indicating a phase transition is detected down 
to 2 K. A Curie-Weiss fitting using MjH = CUT - ©cw) to the data of 0.01 T below 20 K 
gives C = 0.620(1) emu/K mol and ©cw = 0.36(1) K. The small value of ©cw indicates the 
existence of weak ferromagnetic interactions among Ce ions when M H a. 
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Fig. 4. (Color online) (a) Temperature dependences of DC susceptibility for single crystal of CeOBiS 2 mea¬ 
sured at 0.01 and 7 T in // // a. The solid line represents the Curie Weiss fitting line below 20 K. (b) Reciprocal 
susceptibility in // // a (solid circles) and H H c (solid triangles) measured at 1 T and (c) low-temperature mag¬ 
netization in H H a and H H c measured at 2 K. The determined CEF level scheme is shown in (b). The broken 
lines are the calculated results of CEF analysis as described in the text. 
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Table II. CEF energy levels and the corresponding wave functions for CeOBiS 2 . The mixing parameter rj 
varies as \t]\ < 0.13 depending on 


E„(K) 

symmetry 

wave functions 

-710 

^ 7 

Vl -??2|±5> + 77 

+ 5> 

-310 

p(l) 

^ 7 

77|+§)- 


0 

Efi 




Figure 4(b) shows the temperature dependence of the reciprocal susceptibility H/M mea¬ 
sured in 1 T for H H a and H // c. The data demonstrate that Ce 4/-electron magnetic mo¬ 
ments have significant magnetic anisotropy. A shoulderlike anomaly appearing at approxi¬ 
mately 100 K for // // c is a characteristic feature caused by CEF effects. This fact indicates 
clearly that 4/-electrons are in a localized state and the valence of the Ce ions is 3-1-.^^ Fig¬ 
ure 4(c) shows the magnetic field dependence of magnetization at 2 K. The larger values of 
M for H H a compared with those for H H c indicate that the ah-plane is the magnetic easy 
plane. 

Owing to CEF effects with the tetragonal point symmetry C^v (4mm), the J = 5/2 multi- 
plet of a Ce^”^ ion splits into three doublets. The CEF Hamiltonian can be expressed as 
“Kcef = ^ 2^2 ^ 4^4 + where O] are Steven’s operators and B] are CEF parameters. 

A fitting of the CEF model to the H/M(T) data was performed. The best-fit CEF parameters 
are 5^ = 38.1(1) K, 5° = -0.282(4) K, and < 2 K, using which the model calculations ofx 
and M are drawn in Figs. 4(b) and 4(c). Precise determination of (or the mixing parameter 
Tj defined in Table II) is difficult since the behavior of H/M(T) for H 1/ c is insensitive to 
the ^4 term when is small as in the present case. It is remarkable that the shoulderlike 
anomaly appearing in H/M{T) for H 1/ c is well reproduced by the CEF model, confirming 
our interpretation of localized 4/-electrons of Ce^'^ ions. The CEF ground state is a pure = 
+ 1/2 doublet, and two excited doublets are located with the energy separations Ai = ~ 310 K 
and A 2 = ~ 710 K. The wave functions of each level are listed in Table II. 

The 4/-electron contribution to the specific heat C 4 / of CeOBiS 2 was obtained using 
C 4 / = CceOBiS 2 “ C'LaOBiS 2 > whcrc CLaOBiS 2 is the specific heat of EaOBiS 2 single crystals mea¬ 
sured in 0 T. The resulting data are plotted as C 4 //r vs log T in Fig. 5(a). No anomaly indi¬ 
cating a phase transition was detected down to 0.2 K in any field, indicating that 4/-electrons 
are in a paramagnetic state in the measured T and H regions. The most salient feature in 
Fig. 5(a) is that C^f/T in zero field continues to increase with decreasing temperature below 
~ 3 It shows distinct - log T dependence for about one decade in temperature. 
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Fig. 5. (Color online) Temperature dependences of (a) C 4 //r and (b) ^ 4 / for single crystals of CeOBiSa at 0 
and 8 T in // // a. (a) Inset: Temperature dependence of C 4 / at 8 T in // // a. The solid line is the fitting curve 
using the two-level Schottky model. The absolute value of 54 /( 7 ’, 0 T) is determined using the Maxwell relation 
dMldT\H = dSldH\T in combination with 5 4 / (0.3 K, 8 T) = 0 and measured M{T, H) data. 


The 4/-electron contribution to the entropy S 4 / calculated using the C 4 / data is shown in 
Fig. 5(b) as a function of logT. In the temperature range of around 5 - 30 K, S 4 / shows 
plateau behavior with a value of R In 2, consistent with the CEF level scheme obtained above. 
From this fact, it is obvious that the - log T dependence of C^f/T originates in the = +1/2 
ground-state doublet. The characteristic temperature T* of the - log T behavior can be es¬ 
timated from the fitting using C 4 //r oc log(r jT*), which gives T*= 1.8(1) K. Therefore, a 
measurement of the magnetic susceptibility below T* is needed to investigate the magnetic 
properties related to the - log T behavior; in Fig. 4(a), the deviation from the Curie-Weiss 
behavior is hardly noticeable above 2 K. 

By applying fields, the - log T divergence of C^fjT is suppressed and it shifts to higher tern- 
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peratures with its structure changing. Eventually in 8 T, a Schottky peak structure is observed 
at approximately 6 K. This peak is attributable to the thermal excitation in the Zeeman-split 
ground-state doublet, reflecting the fact that quantum fluctuations in 4/ magnetic moments 
dominating in zero field are significantly suppressed in 8 T. The peak height of 3.60 J/K 
mol, which is slightly smaller than 3.65 J/K mol expected for a two-level system, reflects the 
broadening of the peak probably caused by remaining quantum fluctuations. The Zeeman en¬ 
ergy separation estimated from the peak position is 15.1 K, which is slightly larger than 13.8 
K calculated for the CEF ground state; this indicates the existence of weak ferromagnetic 
interactions among Ce ions. 

Despite the nearest-neighbor (nn) Ce-Ce distance being only 3.75 A, no magnetic ordering 
occurs down to 0.2 K. This is attributable to several characteristic features of the crystal 
structure (see Fig. 1). CeO block layers are separated from each other by a large distance of 
c (= 13.49 A), and two BiS 2 layers weakly bound by van der Waals forces are inserted in 
between.This feature probably makes CeO interlayer interactions extremely weak, making 
this material a nearly 2D system (note that, in electron-doped isomorphic Nd(OF)BiS 2 , the 
2D nature of the Fermi surface topology has also been confirmed by Angle-Resolved Photoe¬ 
mission Spectroscopy measurements^^). 

Magnetic interactions in a CeO block layer are probably mediated by superexchange in¬ 
teractions. In a CeO layer, Ce ions form an undulating ichimatsu-patterned (checkerboard- 
patterned) square lattice, which consists of two square sublattices A and B intervening with 
each other. There are two types of Ce-O-Ce bond and each provides different interactions, 
J\ (inter-sublattice) and J2 (intra-sublattice). Because these bonds have the same total bond 
length of 4.69 A and slightly different bond angles of 106 and 116°, and J 2 may com¬ 
pete with each other, causing a CeO layer to become a frustrated magnet.In addition, the 
4/-electron spatial charge distribution of the = +1/2 ground state is elongated along the 
c-axis.^^ Therefore, it is expected that Ji and J 2 each have sufficiently low strength. 

The internal degrees of freedom that the = ±1/2 ground state have are magnetic moments 
of Jx, Jy (Ts), and (r 2 ). Therefore, interactions among Ce ions can be expressed with the 
spin-1/2 XXZ model. In a 2D quantum spin-1/2 XXZ system along with the J\ - J2 frus¬ 
tration, theoretical calculations show that no long-range order (FRO) appears in a certain 
parameter range called the “quantum paramagnetic phase”.Therefore, it may be possi¬ 
ble that CeOBiS 2 is located in this parameter region (or in the vicinity of the QCPs existing 
around it), where quantum fluctuations are expected to dominate at low temperatures close to 
T = 0 . 



J. Phys. Soc. Jpn. 


Considering that ©cw = + 0.36 K reflects the ferromagnetic interaction among Ce ions at the 
wave vector q = 0, the fact that S 4 / starts to decrease below approximately 4 K indicates the 
existence of dominant antiferromagnetic {q 0 ) components with the characteristic energy 
of about 4 K, which causes the - log T divergence of C 4 //r. The lack of magnetic ordering 
down to 0.2 K, which is more than one order of magnitude smaller than 4 K, is in agreement 
with the frustration scenario. 

Among the geometrically frustrated nonmetallic magnets with no LRO, many have been dis¬ 
cussed as candidate of “quantum spin liquids”^^ (e-g-, NiGa 2 S 4 ^^ and Na 4 lr 308 ^'*). These 
materials show roughly -dependent specific heat at low temperatures, and no compounds 
that show exotic - log T dependence except for CeOBiS 2 have been reported in the literature. 
Therefore, the present findings strongly indicate that CeOBiS 2 is the first material found to 
be located at a QCP. 
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